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INTRODUCTION 

The  U.S.  Hydrograph  Laboratory  of  the  Agricultural  Research  Service, 
working  with  the  Central  Technical  Unit  of  the  Soil  Conservation  Service, 
reviewed  a  number  of  existing  hydraulic  and  hydrologic  programs  for  IBM 
1620  digital  computers  to  see  if  any  of  them  might  be  adapted  to  the  water- 
shed project  formulation  needs  of  the  Soil  Conservation  Service  and  hydro- 
logic  research  needs  in  ARS  watersheds.  A  program  developed  by  the  U.S. 
Army  Corps  of  Engineers  [1]  for  computing  water  surface  profiles  by  the 
step  method  [6]  was  selected  and  modified  to  spatially  vary  flow  concor- 
dantly  with  drainage  area,  thus  providing  the  profiles  for  selected  flows- 
per-unit-area . 

Spatial  variations  in  flow  and  drastic  differences  in  cross  sections 
of  the  upland  watershed  streams  used  for  trial  runs  caused  the  computed 
profiles  to  undulate  unrealistically .   A  review  of  results  indicated  that 
undulations  were  greatest  in  reaches  having  steep  slopes.   As  a  result  the 
program  was  further  modified  to  interpolate  hydraulic  elements  for  theoreti- 
cal cross  sections  at  predetermined  vertical  intervals.   Spatial  variation 
of  flow  was  also  graduated  by  interpolating  flow  between  cross  sections  in 
the  ratio  of  the  accumulative  interval  to  the  total  interval  of  elevation. 
By  these  means  the  water  surface  profiles  were  stabilized  in  close  conform- 
ity with  profiles  computed  by  more  sophisticated  techniques  [3], 


COMPUTER  NEEDS 

This  program  is  adapted  for  a  1620  computer,  40K  storage  capacity  with 
a  card  reader,  a  card  punch,  and  an  IBM  407  tabulating  machine  with  an  80-80 
control  panel  for  listing  punched  output.   This  program  is  written  in  Fortran 
II, 

PROGRAM  FEATURES 

The  analytical  procedure  is  basically  an  iterative  solution  of  Ber- 
noulli's Energy  Theorem  [5],   Solutions  are  obtained  to  the  nearest  0.1  foot. 
Friction  losses  are  obtained  from  Manning's  formula.   Bridge  losses,  bend 
losses,  and  contraction-expansion  losses  are  determined  by  the  entry  of  coef- 
ficients [2,7]  at  the  appropriate  cross  sections. 


Output  from  the  program  is  a  listing  of  water  surface  profile  eleva- 
tions for  each  csm  assumed.   By  interrogating  console  sense  switch  No.  3, 
data  for  plotting  storage  indication  curves  will  be  computed  and  summarized 
in  tabular  form. 


BASIC  RELATIONS  FOR  COMPUTATIONS 
Bernoulli ' s  Theorem 

E1    =  E2  +  hv2  +  hf  -  hvx 
wherein: 

Ex  is  elevation  water  surface  upstream  end  of  reach,  in  feet; 
hvx  is  velocity  head  upstream  end  of  reach,  in  feet; 

E2  is  elevation  of  water  surface  downstream  end  of  reach,  in  feet; 
hVg  is  velocity  head  downstream  end  of  reach,  in  feet; 

hf  is  friction  head  loss,  in  feet. 


Manning ' s  Formula 

hf  =  (  ,n?Q  P/a  T  x  RCH 

r   V  1.486  AR2/3  J 

wherein: 

Q  is  discharge  at  section,  in  cubic  feet  per  second; 

A  is  cross-sectional  area,  in  square  feet; 

R  is  hydraulic  radius,  in  feet; 

n  is  Manning's  "n"; 

RCH  is  reach  length,  in  feet. 

Bridge  Head  Loss  (from  Yarnell  [7]  and  Chow  [2]) 

V2 
HBR  =  2K  (K  +  lOuj  -  0.6)  (a  +   15a4)  T" 


2g 


wherein: 


K  is  Yarnell' s  bridge  coefficient  (entered  from  Yarnell); 

to   is  ratio  of  the  velocity  head  immediately  downstream  of  bridge 

to  the  depth  of  flow  through  bridge  (computed  in  program) ; 
<y     is  contraction  ratio  based  on  ratio  of  areas  (computed  in 

program) ; 
V  is  velocity  head  in  feet  (computed  in  program) ; 
g   is  acceleration  due  to  gravity. 


The  values  of  K  for  the  principal  pier  shapes  tested  by  Yarnell  [7] 
are  listed  below. 

Semicircular  nose  and  semicircular  tail  -  -  -  -  0.90 

Twin-circular  piers  with  connecting  diaphragm  -  0.95 

Twin-circular  piers  without  diaphragm  -----  1.05 

90°  triangular  nose  and  90°  triangular  tail  -  -  1.05 

Square  nose  and  square  tail  ----------  1.25 

"Piers  with  one  or  both  ends  lens- shaped  or  convex- shaped  appear 
to  have  a  coefficient  very  nearly  the  same  as  that  for  one  with 
semicircular  ends." 

The  following  contraction  and  expansion  loss  coefficients  [1]  (SQZE) 
are  suggested  for  use  in  the  equation: 


SQZE 


V  2_v  2 
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For  converging  reaches: 

gradual  convergence 
abrupt  convergence 

For  diverging  reaches: 

gradual  expansion 
abrupt  expansion 


0.10 
0.50 


0.20 
0.50 


Bend  loss  coefficients  [4,8,9]  (BEND)  for  use  in  the  equation  hL  = 
BEND  (V12+  V22/  2g)  can  be  computed  as  0.20  (stream  width/radius  of  curva- 
ture) for  all  bends  36°  to  180°. 


FIELD  SURVEY  FOR  BASIC  DATA 

Many  assumptions  and  averages  of  the  hydraulic  characteristics  of  a 
stream  system  are  necessary  to  facilitate  the  computation  of  water  surface 
profiles.   Hence,  it  is  important  that  the  field  survey  document  those  fea- 
tures and  changes  in  features  of  the  stream  and  valley  which  tend  to  accel- 
erate, retard,  contract,  or  expand  the  flow  of  water. 

The  following  should  be  observed  in  obtaining  stream  and  valley  cross 
sections. 


1.  Cross  sections  are  needed  at  points  where  the  stream  flow 

gradient  changes. 

2.  Cross  sections  are  needed  at  points  of  contraction  or  expan- 

sion of  the  channel. 

3.  Cross  sections  are  needed  in  tributaries  immediately  above 


the  confluence  and  in  the  mainstem  immediately  below  the 
confluence . 

4.  Cross  sections  are  needed  to  represent  Manning's  n,  i.e., 

take  sections  so  that  there  will  be  a  uniform  transition 
between  sections,  avoiding  abrupt  changes  in  n  value 
whenever  possible. 

5.  At  bridges  considered  as  controls,  take  cross  sections  at 

center  of  the  bridge  and  immediately  upstream  and  downstream 
of  the  bridge.   (Sufficient  field  notes  should  be  taken  to 
determine  Yarnell's  bridge  loss  coefficient  [2,71). 

6.  Cross  sections  are  required  immediately  above  and  below  rated 

control  sections,  such  as  weirs,  flumes,  etc. 

7.  Distance  between  cross  sections  normally  should  not  exceed 

one  mile.   Shorter  reach  lengths  are  required  in  areas  of 
very  steep  bed  slope.   Determine  both  the  channel  and 
valley  length  between  cross  sections. 

8.  The  horizontal  reference  must  always  be  to  the  left  of  the 

first  cross  section  coordinate  point.  The  survey  must  be 
referenced  consistently,  either  facing  upstream  or  facing 
downstream. 

9.  Cross  sections  must  extend  vertically  to  exceed  the  antici- 

pated elevation  of  the  maximum  flow  profile  to  be  computed. 


DATA  TABULATION 

The  spread  sheet  of  Table  1  is  a  complete  array  of  data  needed  at 
each  cross  section.   The  column  numbers  across  the  top  of  the  table  are  in 
reference  to  the  standard  punch  card  and  the  numbers  listed  in  column  80 
are  the  serial  numbers  of  the  punch  cards  on  which  data  and  computations 
from  the  field  survey  are  to  be  entered.  This  form  is  an  instruction  to 
the  card  punch  operator  but  it  also  serves  as  a  guide  in  the  field  survey. 
For  example,  cards  numbered  3-10  indicate  that  coordinate  points  are  limi- 
ted to  32  for  any  one  cross  section.  There  are  no  limits  concerning  depres- 
sions on  the  flood  plains,  i.e.,  secondary  channels  may  be  topographically 
defined  and  will  be  included  as  part  of  overbank  flow. 

The  lower  portion  of  Table  1  provides  information  and  parameter  cards 
for  eight  additional  profiles  using  the  same  cross  section,  depicted  on 
cards  3-10.   Information  cards  (No.  1)  may  differ  in  reach  length  (columns 
7-11),  and  possibly  in  the  loss  coefficients  (columns  20-28),  for  the  vari- 
ous csm  profiles.   Parameter  cards  (No.  2)  will  differ  in  starting  water 
surface  elevation  (columns  22-26)  at  the  starting  section  only.   The  No.  2 
cards  remain  the  same  for  all  profiles  for  cross  sections  above  the  starting 
section.   However,  they  are  duplicated  for  each  profile  to  avoid  confusion 
in  machine  loading. 


CARD  PUNCHING 

A  list  of  symbols  is  given  at  the  end  of  the  text .  The  card  column, 
format,  and  symbolic  name  of  the  associated  data  fields  are  given  in  Table 
2.   Data  from  Table  1  should  be  punched  as  follows: 

1.  Prepare  only  one  profile  header  card  for  each  water  surface 

profile  to  be  determined. 

(a)  SLNEW  and  HVNEW  columns  will  always  be  blank  except 

when  needed  in  (c) . 

(b)  EI  is  usually  taken  as  1  foot  for  a  first  trial. 

(c)  SLNEW  and  HVNEW  values  are  entered  when  determining 

water  surface  profiles  for  tributaries  and  when  it 
is  desired  to  finish  a  partially  completed  profile. 
Table  4  is  an  example  of  data  prepared  for  a  tribu- 
tary with  its  confluence  between  cross  sections  14 
and  15  of  the  mainstem.   The  starting  section  was 
section  15  (mainstem).   The  SLNEW  and  HVNEW  values 
were  obtained  from  Table  6  (output  data)  in  the 
Appendix. 

2.  Prepare  one  cross  section  information  card  for  each  profile 

at  each  cross  section. 

(a)  If  there  is  no  left  overbank  area,  OBL  must  be  alge- 

braically less  than  the  first  x-distance,  e.g., 
x-distance  100  feet;  OBL  may  be  99  or  less. 

(b)  If  there  is  no  right  overbank  area,  OBR  must  be 

algebraically  greater  than  the  longest  x-distance. 

(c)  Yarnell's  bridge  loss  coefficients,  bend  loss  coef- 

ficients, and  contraction-expansion  loss  coeffi- 
cients are  taken  from  references  [1,4,8,  and  9]. 

(d)  Cross  sections  are  identified  (columns  1-4)  with 

numeric  characters  beginning  with  one  for  the  up- 
stream section. 

(e)  Streams  are  identified  (columns  5-6)  with  two  alpha- 

betic characters  (e.g.  SB  for  South  Branch). 

30  Prepare  one  cross  section  parameter  card  for  each  profile  at  each 
cross  section. 

(a)  When  possible  the  starting  water  surface  elevation  for  a 

given  csm  is  obtained  from  existing  rating  curves. 

(b)  For  resumed  computations  or  for  tributary  profiles, 

given  in  step  1(c),  the  starting  ELEV  becomes  the 
elevation  previously  computed  for  that  water  surface 
profile  at  that  section  (cross  section  15  of  example) . 


TABLE  2. --Input  Data  Card  Format* 


Profile  Header  Card; 

1-3  4-9  10-16  17-20 
X.XX  XXXXXX.  .XXXXXXX  XX. XX 
EI    CSM      SLNEW      HVNEW 

Note:   The  profile  header  card  is  coded  with  a  zero  in  column  80. 


Cross  Section  Information  Card: 

1-4  5-6  7-11  12-15  16-19  20-22  23-25  26-28  29-33  34 
XXXX  XX  XXXXX.  XXXX.  XXXX.  X.XX  X.XX  X.XX  XXXX.X  X. 
NSTA  LET  RCH     OBL    OBR    YARK   BEND    SQZE   EVRT    FLAG 

Note:   The  cross  section  information  card  is  coded  with  a  1  in  col- 
umn 80 . 


Cross  Section  Parameter  Card : 

1-4  5-6  7-12  13-15  16-18  19-21  22-27  28-31 
XXXX  XX  XXXX. XX  .XXX  .XXX  .XXX  XXXX. XX  XX. XX 
NSTA  LET  DA       C(l)    C(2)    C(3)   ELEV     DT 

Note:   The  cross  section  parameter  card  is  coded  with  a  2  in  column 
80. 


Valley  Cross  Section  Coordinate  Data: 

1-4   5-6  7-11    12-16   17-21   22-26   27-31   32-36   37-41 
XXXX  XX   XXXX.X  XXXX.X  XXXX.X  XXXX.X  XXXX.X  XXXX.X  XXXX.X 
NSTA  LET  X(I)    Y(I)     X(I+1)   Y(I+1)   X(I+2)   Y(I+2)   X(I+3) 

42-46  47 
XXXX.X  X. 
Y(I+3)   EX 

Note:   Valley  cross  section  coordinate  data  cards  (Maximum  of  8)  are 
coded  from  3-10  in  column  80. 


*Symbols  are  taken  from  list  of  program  symbols  given  at  the  end  of  the  text 


(c)  DT  must  always  be  entered  for  the  first  downstream  cross 
section  if  storage  indication  data  are  desired.   DT 
entries  are  required  for  subsequent  sections  only  if 
routing  time  is  to  be  varied. 

Only  one  set  of  valley  cross  section  coordinate  data  is  needed, 
Table  3  in  the  Appendix  illustrates  the  data  required  for 
determining  the  10  csm  water  surface  profile  of  Fall  River, 
Kansas.  To  determine  other  profiles  replace  card  types  0, 
1,  and  2;  these  are  the  substitute  cards  of  Table  1.  An 
abbreviated  example  of  input  data  for  a  tributary  profile 
is  given  in  Table  4  of  the  Appendix.   Table  5  illustrates 
input  data  for  a  bridge  section. 


OPERATING  INSTRUCTIONS 

1.  Load  object  deck  with  library  floating  point  subroutines 

(Fortran  II) . 

2.  Turn  switch  3  on  for  storage  indication  data. 

3.  Check  to  see  that  order  of  input  corresponds  with  that  of 

Table  1  or  Table  3. 

4.  For  the  first  section  all  three  control  cards  (card  types  0, 

1,  2)  must  be  present  in  the  order  0,  1,  2. 

5.  For  subsequent  sections  in  the  same  profile  run,  only  two 

control  cards  are  needed  (1,2).  A  card  type  2  is  not  re- 
quired when  a  zero  appears  in  column  34  of  card  type  1. 

6.  When  stacking  a  new  profile  run  key-punch  a  new  profile 

header  card. 

7.  Substitute  card  types  1,  2  with  the  new  control  cards  for 

the  next  csm. 

8.  At  the  completion  of  a  profile  sort  storage  indication  data 

from  water  surface  profile  data.   Storage  indication  data 
are  coded  with  a  _1  in  column  80. 

9.  Exceeds  survey:  When  the  survey  bounds  are  exceeded  the  instruc- 

tion EXCEEDS  SURVEY  LIMITS  (section  at  which  boundaries  were 
exceeded)  TURN  SW  1  ON  FOR  RERUN  WITH  SHORTER  REACH  LENGTH 
OFF  CONTINUES  UPSTREAM  is  typed  by  the  typewriter.   Check  the 
EI  interval.   Turn  SW  1  on  and  push  start.  Enter  one-half 
the  previous  EI  value  by  typewriter  with  format  X.XX.   Turn 
SW  1  off  to  revert  to  automatic  selection  of  EI.   If  the 
water  surface  again  exceeds  survey  for  the  cross  section, 
leave  SW  1  off  and  push  start  .to  continue  with  the  next 
cross  section. 
10.  Partially  completed  profile:  When  a  profile  has  not  been  com- 
pleted due  to  machine  failure  or  other,  the  profile  may  be 
continued  at  the  computed  section  by  key  punching  new  card 
types  0,  1,  2  for  the  last  cross  section  at  which  the  water 
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surface  elevation  has  been  determined.   Enter  the  same  data 
as  required  for  tributary  profiles  (Table  4) .   The  new 
values  to  be  punched  are:   SLNEW,  HVNEW  (card  type  0),  RCH=0 
(card  type  1) ,  and  ELEV,  DT  (card  type  2) . 

OUTPUT  DATA 

An  example  of  output  data  is  given  in  Table  6.   This  output  resulted 
from  the  input  data  of  Table  3.   Flow  through  each  section  is  given  in  cfs; 
starting  water  surface  elevation  is  given  in  feet  and  Manning's  n  values 
are  for  left  overbank,  channel,  and  right  overbank  areas,  respectively. 

Starting  at  cross  section  22W  with  reach  length  (RCH=0)  the  data, 
from  left  to  right,  are  coded  as  follows: 

Example  (Table  6) 

22W 
0.0007845 
0.00 
3.42 
0.00 
0.18 
989.80 

0.00 

0. 
926. 

0. 

0.00 
989.80 

0.00 
7.05 
0.00 
0.00 

It  is  suggested  that  a  clear  plastic  overlay  be  made  for  easy  iden- 
tification of  water  surface  profile  output  data. 

An  example  of  storage  indication  data  is  given  in  Table  7  of  the  Ap- 
pendix.  These  data  were  also  derived  from  the  input  data  given  in  Table  3 
of  the  Appendix.   The  data  appearing  in  the  columns  headed  by  S  and  (0/2  + 
S/DELT)  are  given  as  a  power  of  10,  e.g.  the  first  value  under  S  is 
.38522213  E  +  04,  which  is  equal  to  3852.2213.   If  the  exponent  +04  had 
been  -04,  then  the  decimal  would  be  moved  four  places  to  the  left  or 
.00003852213. 


1st  line: 

Code 

1. 

NSTA,  LET 

2. 

SLOPE 

3. 

V(l) 

4. 

V(2) 

5. 

V(3) 

6. 

HV1 

7. 

ECOMP 

2nd  line: 

1. 

HFRIC 

2. 

A(l) 

3. 

A(2) 

4. 

A(3) 

5. 

BENO 

6. 

EWS 

3rd  line: 

1. 

R(D 

2. 

R(2) 

3. 

R(3) 

4. 

SQZO 
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The  area  given  in  the  last  data  column  is  the  cross-sectional  area  of 
flow.   It  is  listed  for  use  with  flow  rate  to  derive  the  average  velocity 
of  flow,  sometimes  required  for  deriving  travel  time  or  coefficients  for 
certain  types  of  flood  routing. 


Symbol 

Input : 

BEND 

C(I) 
1=1,3 

CSM 

DA 
DT 
EI 

EINT 
ELEV 
EVRT 
EX 

FLAG 

HVNEW 

LET 

NSTA 


SYMBOLS  USED  IN  THE  PROGRAM 

Description 


Bend  loss  coefficient. 

Manning's  n  for  left  overbank,  channel,  and  right  overbank 
areas,  respectively. 

The  flow  rate  in  csm  at  which  water  surface  profile  is  to 
be  determined. 

Drainage  area  above  cross  section  in  square  miles. 

Routing  time  increment  in  hours. 

Elevation  increment.   This  value  determines  the  number  of 
subdivisions  between  cross  sections. 

EI  (may  be  entered  by  typewriter) . 

Starting  water  surface  elevation. 

Invert  elevation  immediately  downstream  of  bridge. 

One  digit  integer  having  the  value  of  0  or  1.  0  if  more 
coordinate  data  for  given  cross  section  to  be  read.  1 
if  last  card  of  coordinate  data. 

One  digit  integer  having  the  value  of  0  or  1.  0  if  no  cross 
section  parameter  to  follow.  1  if  cross  section  parameter 
card  to  follow. 

Velocity  head  of  cross  section  from  which  tributary  pro- 
file computations  begin. 

Stream  or  river  identification.   Two  characters  of  the 
alphabet  may  be  used. 

Cross  section  identification  -  all  cross  sections  must 
have  a  numeric  designation  beginning  with  one  for  the 
upstream  section. 
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SYMBOLS  USED  IN  THE  PROGRAM  (Continued) 

Symbol 

OBL 


Description 


Distance  from  horizontal  reference  point  to  left-most  edge 
of  stream  channel. 


OBR  Distance  from  horizontal  reference  point  to  right-most  edge 
of  stream  channel. 

RCH  Reach  length  -  distance  between  cross  sections. 

SLNEW  Water  surface  slope  of  the  cross  section  from  which  tribu- 
tary profile  begins. 

SQZE  Contraction  -  expansion  loss  coefficient. 

X(I)  x-distance  (horizontal  distance  from  reference  to  the  coor- 

1=1,32  dinate  points  in  the  cross  section). 

Y(I)  Elevation  of  coordinate  points  (normally  referenced  to  mean 

1=1,32  sea  level) . 

YARK  Yarnell's  bridge  loss  coefficient. 

Output: 

A(I)  Area  of  left  overbank,  channel,  and  right  overbank, 

1=1,3  respectively. 

Al  Area  of  cross  section  being  processed. 

A2  Area  of  last  cross  section  analyzed. 

BENO  Bend  head  loss. 

C(I)  See  input. 
1=1,3 

DELT  Routing  time  interval  (hrs.)  -  DT 

ECOMP  Computed  water  surface  elevation. 

ELEV  See  input. 

EWS  Estimated  water  surface  elevation. 

FLOW  Flow  rate  (cf s)  . 

FL0W1  Flow  rate  (cfs)  at  last  cross  section  analyzed. 
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SYMBOLS  USED  IN  THE  PROGRAM  (Continued) 

Symbol  Description 

FL0W2  Flow  rate  (cfs)  at  cross  section  being  processed. 

HBR  Bridge  head  loss. 

HFRIC  Friction  head  loss. 

HV1  Velocity  head  at  section  being  analyzed. 

KODE  Code  1  punched  in  column  80. 

LET  See  input. 

NCSM  Flow  rate  in  csm. 

NSTA  See  input. 

NSTAO  Cross  section  identification  of  last  section  analyzed. 

R(I)  Hydraulic  radius  of  left  overbank,  channel,  and  right  over- 
1=1,3         bank,  respectively. 

RCfi  See  input. 

R0F1  FL0W1/2  +  S/DELT 

S  Storage  between  cross  sections  (cfs-hrs.). 

SLOPE  Water  surface  slope. 

SQZO  Expansion-contraction  head  loss. 

V(I)  Average  velocity  (ft/sec)  of  left  overbank,  channel,  and 
1=1,3         right  overbank  respectively. 


Main  Program; 


AA(I) 
1=1,3 

ALFA 

AL4 

AN 

AREA 


Area  of  left  overbank,  channel,  and  right  overbank  of  last 
section  analyzed. 

Channel  constriction  ratio. 

ALFA  raised  to  the  4th  power. 

A  count  of  the  number  of  intermediate  cross  sections  processed, 

Sub-program  for  determining  the  area  of  cross  sections. 
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SYMBOLS  USED  IN  THE  PROGRAM  (Continued) 

Symbol 

ASUM 


Description 


Sum  of  the  left  overbank,  channel,  and  right  overbank  areas, 
respectively. 

AVGA        Average  of  two  end  areas. 

CC(I)       Manning's  n  of  the  left  overbank,  channel,  and  right  overbank 
1=1,3         for  the  last  cross  section  analyzed. 

CFAC  Number  of  reach  sub-divisions  between  cross  sections. 

CHE  Elevation  difference  between  cross  sections. 

CHB  Elevation  of  channel  bottom  at  cross  section  being  processed. 

CHBO  Elevation  of  channel  bottom  at  last  cross  section  processed. 

DELHV       Difference  in  velocity  heads  between  two  cross  sections. 

DELQ(I)     0.01486  A(I)  R(I) ^  ^   discharges  based  on  sk   =  0#01# 

1=1,3  nK) 

ECOMO  Elevation  of  water  surface  at  -last  cross  section  analyzed. 

ECOMO0  ECOMO 

ETEST  Test  for  balance  of  energies. 

EWS2  Estimated  water  surface  elevation  for  intermediate  sections. 

HVO  Velocity  head  at  last  section  processed. 

I  Index  of  do  loop. 

II  1+3  (determines  upper  limit  of  read  loop) . 

J  Index  of  do  loop. 

JJ  Index  code  having  the  value  of  0  or  1.   0  if  channel  elevation 

to  be  determined.   1  if  channel  bottom  elevation  for  given 
section  has  been  determined. 

K  Index  code  having  the  value  of  1  or  2.   1  for  first  cross 

section  and  bridge  computations.   2  for  remainder  of  given 
profile  analyses. 
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SYMBOLS  USED  IN  THE  PROGRAM  (Continued) 

Symbol 

LL 


Description 


LN 

LN1 

M 
N 
NPASS 

NSECS 
OMEGA 

P(D 
1=1,3 

QSUM 

RATIO 

RCHO 

RR(D 
1=1,3 

SCH 

SLOPO 

SLOPOO 

SUMVQ 


Index  code  having  the  value  of  0,  1,  or  2.   0  to  continue 
on  normal  flow  path.   1  for  shorter  reach  length,  then 
branch  back  to  beginning  of  section  computations.   2  to 
proceed  to  next  section  (beyond  survey  bounds) . 

Index  code  having  the  value  of  0,  1.   0  if  elevation  dif- 
ference is  not  to  be  determined.   1  if  elevation  differ- 
ence to  be  determined. 

Index  code  having  the  value  of  0  or  1.   0  for  first  cross 
section.   1  for  subsequent  sections. 

1+3  (determines  upper  limit  of  search  loop)  . 

Counts  number  of  coordinates  per  cross  section. 

Index  code  -  counts  number  of  times  estimated  water  surface 
exceeds  survey  bounds. 

CFAC 

Velocity  head  immediately  downstream  of  bridge  divided  by 
depth  of  water  immediately  downstream  of  bridge. 

Wetted  perimeter  of  left  overbank,  channel,  and  right 
overbank,  respectively. 

Sum  of  DELQ(I) . 

Ratio  of  discharge  for  which  water  surface  profile  is  being 
computed  to  a  discharge  based  on  S2  =  0.01. 

Distance  between  two  given  cross  sections. 

Hydraulic  radius  of  last  computed  cross  section. 


Channel  slope. 

Water  surface  slope  at  last  computed  section. 

SLOPO  (before  intermediate  sections) . 

Sum  of  products  of  velocities  squared  times  discharges  for  left 
overbank,  channel,  and  right  overbank,  respectively. 
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SYMBOLS  USED  IN  THE  PROGRAM  (Continued) 
Symbol 


Description 


V2Q(I)      V2Q  for  Left  overbank,  channel,  and  right  overbank,  respec- 
1=1,3         tively  -  used  to  determine  weighted  average  of  stream 
velocity. 

YDELO       Elevation  of  preceding  coordinate  point. 

Sub-program; 


A  (I) 

1=1,3 

AN 

AREA 

B 

C 

D,E,G,H 

EVRT 

EWS 

EWS2 

I 


LET 

LL 

LM 

N 
NPASS 


See  output. 

See  main  program. 

See  main  program. 

Horizontal  distance  between  coordinate  points. 

Vertical  distance  between  coordinate  points. 

Arguments  of  arithmetic  function  statements. 

See  input . 

See  output. 

See  main  program. 

Index  of  do  loop. 

Index  code  having  value  of  1,  2,  or  3.   1  for  left  overbank. 
2  for  main  channel,  and  3  for  right  overbank. 

Index  code  for  routing  to  correct  arithmetic  function 
statement. 

See  input . 

See  main  program. 

Index  code  for  determining  limits  of  EWS2  when  estimate  ex- 
ceeds survey  bounds . 

See  main  program. 

See  main  program. 
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Symbol 

NSTA 

See 

input . 

OBL 

See 

input . 

OBR 

See 

input . 

P(I) 

See 

main  program. 

1=1,3 

RCH 

See 

input . 

RCHO 

See 

main  program. 

SCH 

See 

main  program. 

SLOPOO 

See 

main  program. 

X(I) 

See 

input . 

1=1,32 

xws 

Y(I) 
1=1,32 

YDEL 

YDELO 


Description 


Distance  from  horizontal  reference  point  to  edge  of  water 
surface. 

See  input . 


Depth  of  water  at  coordinate  point  being  processed. 
See  main  program. 
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APPENDIX 


TABLE   3.   INPUT  DATA  FOR  DETERMINING  THE  WATER  SURFACE  PROFILE 
OF  THE  MAIN  STEM  OF  UPPER  FALL  RIVER.  KANSAS  AT  A  FLOW  OF  10  CSM. 

100000010  0 

0022  W0000030163185000000000000001  10        1 

0022  W0317380620450500989800025  10         2 

0022  W000001006301«+5010024025101000705350099960  3 

0022  W07930099661185009991138500999416200099900  4 

0022  W20050099872975009986301000998530160099360  5 

0022  W30350098873045009813315000982031600099700  6 

0022  W31850100253220010027345501004840400100641  7 

0021  W1510032003520000000000000001  10        1 

0021  W03107306204507000000  10         2 

0021  W00000103580050010325022501022304050101810  3 

0021  W07900101690920010168115001016311600101600  4 

0021  W11850101211325010123149501014716650101830  5 

0021  W22700101912510010199277001020431300102130  6 

0021  W32000101933280010103333501008833500100350  7 

0021  W33800099403475009943348501011435150101680  8 

0021  W35200102103865010229409501023443550102510  9 

0021  W4555010273.5120010296521001033253100105091  10 

0020  W1050054805680000000000000001  10        1 

0020  W03042405805005000000  10         2 

0020  W00600104690540010426064001041208200103650  3 

0020  W09650103541555010347179501032920550103260  4 

0020  W24150103102930010288303501028333500103130  5 

0020  W38550103034285010271483501027449550102750  6 

0020  W50250102885155010264530501023753700102310  7 

0020  W54700101915480010176550001013055200100O50  8 

0020  W56100100105615010110562501014256400101900  9 

0020  W5680010469                                  1  10 

0019  W0520050355209000000000000001  10        1 

0019  W02986806304506000000  10        2 

0019  W00400105520390010509073001047108600104380  3 

0019  W1370010401143001040018400103792U50103770  4 

0019  W21700103632325010330267001034827550103880  5 

0019  W36100103623755010333427501034544300102970  6 

0019  W46400102764735010254483501025349350102370  7 

0019  W50350102635065010182509001010451600101030  8 

0019  W51700101745185010226520001034852090103600  9 

0019  W551001037657500104195915010560           1  10 

0018  W0800004450675000000000000001  10        1 

0018  W02118406005006000000  10        2 

0018  W00000105650220010507040501049204250104630  3 

0018  W04450104270475010359048001026506200102610  4 

0018  W06300103600675010410079001037511000104040  5 

0018  W13700104451490010480156001048217550104720  6 

0018  W216001047923150104842365010584           1  7 

0017  W06800-1100354000000000000001  10        1 

0017  W02092800005505000000  10        2 

0017  W00000106980015010671008501045701800104500  3 

0017  W02150104430230010402033501040103450105100  4 

0017  W03540105540360010550048001054105550105580  5 

0017  W07200105740721010576092001048515800104890  6 

0017  W17400105002080010564235001056525400105550  7 

0017  W30750105543250010539338001054335650105480  8 

0017  W373001057239000106154040010707           1  9 

0016  W1330024602660000000000000001  10        1 

0016  W02050005006000000000  10        2 

00  1 6 >W00000 10 84300 250 10 7790050010 7360 1300 106850  3 


20 


TABLE  3.   (CONTINUED) 


0016  WO 380010696 1195010 6981395010 68 719250 
00  16  W2 165 01 064724  300 1062524600 106 18  248  50 
0C16  W2505 01 043 72 5950 1044 126000 10 5 142655  0 
0015  W080001 105 12  5000000000  0000001 
0015  W01930307907010000000 

0015  W008 00 10865 01 350 108220 175 010 789 02050 
0015  W044 5010755058 50 107540  72001071608250 
0015  W0915010 696 09850106731 02 0010 69511050 
00  15  Wl 125 010687 11 300 10658 11450 10572 12250 
00  15  W 1230 01 062 8 12500 10 748 12900 10 738 14250 
0015  W14700107171 570010839160001087000000 
0014  W 10400 182 520 1000000000000000 12 00 
0014  W01097905006000000000 

0014  WOOOOO 11 14101400 11 04702000 11 04503000 
0014  W04600 1099705  700 109700  71 5010963  078  50 
0014  W08 100 1094608 500 1096 10925 01 0975 11 100 
0014  Wl 2900 10963 13850 10949 14650109341 5500 
0014  W17150109031745010 89 51 78 001087718250 
0014  W18350108311857010 75 3192201075519270 
00  14  W 19420 108 56 196 20 10 990 2 00 70 11 170 
0013  W0830006 150 72 000000000000000 12 00 
0013  W01040307507006500000 

0013  WOOOOO 1121 2013501 11600285011 15706150 
0013  W0635 01097606400 10969065  5  0 10 877  07 000 
0013  W071 00 1088607200 109720835 01 100708  700 
0013  W 15000 10978 15 500 10 979 17500 11 0462 0900 
0013  W23 100 1098725400 11 020327001 108 136000 
0012  W0470 01 3 05 142400000000000000 1200 
0012  W00979908306504500000 

0012  W000001126000 5001123 90 250011 17204150 
0012  W06950111321000011125104501107711500 
0012  W1230011 04 113050110261320010 99713220 
0012  W 13820 10929 14 100 11 088 14240 11 139 17000 
0012  W26200 11 14927700 11 15632600 11 162 39700 
00  12  W40500 11 182429001 1222446001 125347050 
0011  Wl  1300-1 100  25600000000000000 12  00 
0011  W00946400007005000000 

0011  W011 00 11 379014001 12660 1550 11 19701  700 
0011  WO  1740 11 07702400 11 08 10 2520 11 203 0  2 560 
0011  W0295 01 123 103700 11 27305850 11 28507 100 
0011  W086501 1306107001 131 5 14200 11 338 17300 
0010  W09300370 54 16000000000000000 1200 
0010  W00878106007200000000 

0010  WO  1300 1169702950 11 66705200 11 63 008400 
0010  W09900 11600 14 100 11 58 31 8 100 11 5552 1800 
0010  W238001 15242840011 51 13070011 503  32600 
0010  W34900 11 46935 100 11 4683 5200 11 45937050 
0010  W37450 11 370 375701 135 138150 11 35038250 
0010  W384001 1459388001 152939900 1160840650 
0010  W4150011729 

0009  W05900 18 10205000000000000000 1200 
0009  W00572107108500000000 

0009  WOOOOO 1175702 0001 171 50 5 700 11648 08500 
00  0  9  W107001 1580 114001 157 11 3250 11 55 11 3700 
0009  W 15 100 11 580 15250 11 568 16300 11 538 17000 
0009  W1715011 52 018100 115331 8200 1151918250 
00  09  W18780 1143818950 11 51 3 19050 11 57620450 


106810 
105430 
108971 


107580 
107150 
107110 
105770 
107360 
000001 


110310 
109510 
109760 
109050 
108520 
107870 
1 


111340 
108780 
110370 
110230 
112121 


111340 
110420 
109220 
111410 
111660 
112951 


111580 
112340 
112890 
113981 


116100 
115340 
114870 
114650 
114360 
116610 
1 


116080 
115660 
115370 
114420 
118581 


10 
10 


10 
10 


10 
10 


10 
10 


10 
10 


10 
10 


10 
10 
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TABLE  3.   (CONTINUED) 

0008  W0600009951218000000000000001200  10        1 

0008  W00536810508515000000  10        2 

0008  W00600118800145011802017501178303150117600  3 

0008  W04250117490620011716069501167807650116870  4 

0008  W08050116930855011667090501166609350116600  5 

0008  W09950116591015011634109501161311350116100  6 

0008  W11750116091200011591121301159012180117340  7 

0008  W12240117331244011734127401176014090117850  8 

0008  W1429011880                                  1  9 

0007  W0370010251140000000000000001200  10         1 

0007  W00388604007015000000  10        2 

0007  W00000121140150012064020001198103400119310  3 

0007  W04150119010171011881091001187109800118320  4 

0007  W10100118111025011797107001174911050117100  5 

0007  W11350117111140011762121001176412550117650  6 

0007  W12650118181270011889129501199613100121141  7 

0006  W03700-1100154000000000000001200  10        1 

0006  W00361800006504000000  10         2 

0006  W00000121160055011980007001191000850118220  3 

0006  W01500118210154011913017501188802200119060  4 

0006  W02600119400470011944071001194510100119570  5 

0006  W104001201310650120621090012116           1  6 

0005  W0690004100530000000000000001100  10        1 

0005  W00349708205507000000  10         2 

0005  W00000123200030012309008001228201400122370  3 

0005  W0230012215029O01221303250122O00380012168O  4 

0005  W03950121780410012183043301211904630120900  5 

0005  W05050120910511012098051701219805300122210  6 

0005  W05500122260630012227070501224408800122650  7 

0005  W108501229511350123111205012333            1  8 

0004  W0510003950535000000000000001100  10        1 

0004  W00245307505007500000  10'       2 

0004  W02750124850310012473034501244203700123660  3 

0004  W03950123460405012285046501228604950123070  4 

0004  W05150123330535012356058001236608350123670  5 

0004  W10500123851145012400121501243012750124650  6 

0004  W1350012498                                  1  7 

0003  W03900-1100280000000000000001100  10        1 

0003  W00232900006004000000  10        2 

0003  W00500126350095012561014001252601830124860  3 

0003  W01850124510243012450026001248602700125520  4 

0003  W02800125610325012568048001256306500125610  5 

0003  W08500125601040012553139001253815300125370  6 

0003  W16900125681755012635                      1  7 

0002  W06500-1100165000000000000001100  10        1 

0002  W00138700006006500000  10        2 

0002  W00350130300065012853007801274501180127440  3 

0002  W01250127780145012792016501282202150128230  4 

0002  W02450128390450012843060001285107200128540  5 

0002  W08000128671050012874108001288911000130301  6 

0001  W0800008651020000000000000001100  10        1 

0001  W00103805005500000000  10        2 

0001  W00000133240055013272037001321305950131890  3 

0001  W06300131810865013175088001315608900131220  4 

0001  W09150131230920013160093001318809400132540  5 

0001  W09600132921015013324                      1  6 
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TABLE   4  .   INPUT  DATA  FOR  DETERMINING  THE  WATER  SURFACE  PROFILE 
OF  A  TRIBUTARY  TO  THE  MAIN  STEM  OF  UPPER  FALL  RIVER.  KANSAS 

AT  A  FLOW  OF  10  CSM. 


10000001000108100016  10        0 

0015  W0800011051250000000000000001  10        1 

0015  W01930307907010000000  10        2 

0015  W00800108650135010822017501078902050107580  3 

0015  W04450107550585010754072001071608250107150  4 

0015  W09150106960985010673102001069511050107110  5 

0015  W11250106871130010658114501057212250105770  6 

0015  W12300106281250010748129001073814250107360  7 

0015  W14700107171570010839160001087000000000001  8 

0013EB0950006350785000000000000001  10        1 

0013EB0077800500450450000000000001  10        2 

00  13EB00000 11 04800500 1098600700 10962 01 350 109 330  3 

00  13EB022 50 1092 103300 1090 103950 1091306 350 109280  4 
001 3 FB06 55 0108 2 306700 10 76 30675 01 073707500 10735  0  5 
001 3 EB076 00 10792 07850 10 8280 8 300 10 84209200 108 590  6 

001  3 EB 102 50 10863 11000109001 35 5 01 090 11 5400109090  7 
001 3 EB 16600 1092 2 18650 10928 2 225 01093 124 100 109380  8 
00 13EB2 79 5010946 28850109533010010 99730850110130  9 
0013EB32000110263400011048  1  10 
0C12EB10000-1100363000000000000001  10  1 
0012EB0071900000500500000000000001  10  2 
00 12EB000001 12380 18001 1129020001 108502 300109490  3 
001 2 EB02400 108940 3 000 10 8920 32 00 10 9690 36 30 1099 10  4 
00  12EB04300 1099505550 11 00607350 11 0441 0200 110220  5 
00  12EB1 11 001 1028 11700 11 05 11 2000 11 060 13650 110650  6 
00 12EB 1375 01 1058 169501 101 3 17950 110 10 18 050 110 140  7 
00 12EB19 1001 105 120900 1107523250 11 099241501 11020  8 
00  12EB27750 11 103346001 11073625011 12 139450 111500  9 
0012EB4400011238                                1  10 
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TABLE 


INPUT  DATA  FOR  DETERMINING  THE  BACKWATER  EFFECT  OF  A  BRIDGE, 
DATA  FROM  UPPER  FALL  RIVER.  KANSAS 


1000 
0022 
0022 
0022 
0022 
0022 
0022 
0022 
0022 
•0022 
0022 
0022 
0022 
0022 
0022 
0022 
0021 
0021 
0021 
0021 
0021 
0021 
0021 
0021 
0021 
0021 


00010 
W0O0 
W031 

wooo 

W079 

W200 

W303 

W318 

BWOOO 

BW031 

BWOOO 

BW162 

BW190 

BW420 

BW448 

BW750 

W151 

W031 

WOOO 

W079 

W118 

W227 

W320 

W338 

W352 

W455 


0030 
7380 
0010 
3009 
5009 
5009 
5010 
0041 
4050 
0010 
0010 
0010 
5009 
0009 
0010 
0032 
0730 
0010 
0010 
5010 
0010 
0010 
0009 
0010 
5010 


16318500 
62045050 
06301450 
96611850 
98729750 
88730450 
02532200 
20460009 
62045010 
32505800 
02016600 
20025000 
90542800 
95546000 
21078000 
00352000 
62045070 
35800500 
16909200 
12113250 
19125100 
19332800 
94034750 
21038650 
27351200 


0000000000001 

0989800025 

100240  251010007053  50  099  96  0 

09  991138500999416200  099900 

09986301000998530160  0  99  360 

09813315000982031600  099700 

10027345501004840400100  641 

0000000099051 

10295091001021515800102000 

09985168000998317200101200 

10170310001017541200102600 

099104  305009  93044100  09  9400 

10257556001020570100102040 

102258240010330  1 

0000000000001 

00000 

103250  225010223  04  050101810 

10168115001016311600101600 

10123149  501014716650101830 

10199277001020431300102130 

10103333501008833500100350 

09  943348  501011435150101680 

102294  09  501023  44  3550102510 

10296  521001033253100105091 


10 
10 
10 


10 
10 


10 
10 


0 

1 

2 
3 
4 
5 
6 
7 
1 
2 
3 
4 
5 
6 
7 
8 
1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
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TABLE  6    WATER  SURFACE  PROFILE  MACHINE  OUTPUT 
FOR  10  CSM  FLOW  UPPER  FALL  RIVER.  KANSAS 


10  CSM  PROFILE 


FLOW     3173. 

STARTING  W.S. 

ELEV. 

989.80 

N  VALUES     .062 

.045 

.050 

RCH  = 

0. 

22 

W      .0007845 

0.00 

3.42 

0.00 

.18 

989.80 

0.00 

0. 

926. 

0. 

0.00 

989.80 

0.00 

7.05 

0.00 

0.00 

FLOW     3107. 

N  VALUES     .062 

.045 

.070 

RCH  = 

15100. 

21 

W      .0008551 

0.00 

3.48 

0.00 

.  18 

1002.22 

1.07 

0. 

891. 

0. 

0.00 

1002.22 

0.00 

6.78 

0.00 

0.00 

FLOW     3042. 

N  VALUES     .058 

.050 

.050 

RCH  = 

10500. 

2C 

W      .0006001 

0.00 

3.23 

0.00 

.16 

1010.11 

1.07 

0. 

940. 

0. 

0.00 

1010.  13 

0.00 

7.93 

0.00 

0.00 

FLOW     2986. 

N  VALUES     .063 

.045 

.060 

RCH  = 

5200. 

19 

W      .0025314 

0.00 

5.20 

0.00 

.42 

1017.07 

1.44 

0. 

573. 

0. 

0.00 

1017.06 

0.00 

5.51 

0.00 

0.00 

FLOW     2118. 

N  VALUES     .060 

.050 

.060 

RCH  = 

8000. 

18 

W      .0018665 

0.00 

3.56 

0.00 

.19 

1030.45 

.98 

0. 

593. 

0. 

0.00 

1030.44 

0.00 

3.94 

0.00 

0.00 

FLOW     2092. 

N  VALUES    0.000 

.055 

.050 

RCH  = 

6800. 

17 

W      .0022843 

0.00 

3.45 

0.00 

.18 

1045.21 

1.09 

0. 

606. 

0. 

0.00 

1045.18 

0.00 

3.22 

0.00 

0.00 

FLOW     2050. 

N  VALUES     .050 

.060  0 

.000 

RCH  = 

13300. 

16 

W      .0001375 

0.00 

1.70 

0.00 

.04 

1055.26 

.87 

0. 

1204. 

0. 

0.00 

1055.31 

0.00 

9.13 

0.00 

0.00 

FLOW     1930. 

N  VALUES     .079 

.070 

.100 

RCH  = 

8000. 

15 

W      .0025612 

0.00 

4.51 

0.00 

.31 

1062.32 

1.63 

0. 

427. 

0. 

0.00 

1062.37 

0.00 

4.41 

0.00 

0.00 

FLOW     1097. 

N  VALUES     .050 

.060  0 

.000 

RCH  = 

10400. 

14 

W       .0016101 

0.00 

3.26 

o.oo 

.16 

1079.89 

.92 

0. 

336. 

0. 

0.00 

1079.88 

0.00 

3.84 

0.00 

0.00 
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TABLE   6      (CONTINUED) 

FLOW     1040. 

N  VALUES     .075 

.070 

.065 

RCH  = 

8300. 

13 

W      .0014826 

0.00 

3.37 

0.00 

.17 

1092.84 

1.02 

0. 

307. 

0. 

0.00 

1092.84 

0.00 

4.31 

0.00 

0.00 

FLOW      979. 

N  VALUES     .083 

.065 

.045 

RCH  = 

4700. 

12 

W      .0007784 

0.00 

2.62 

0.00 

.10 

1098.25 

.93 

0. 

373. 

0. 

0.00 

1098.27 

0.00 

4.84 

0.00 

0.00 

FLOW      946. 

N  VALUES    0.000 

.070 

.050 

RCH  = 

11300. 

11 

W      .0014732 

0.00 

3.13 

0.00 

.15 

1112.26 

1.  11 

0. 

302. 

o. 

0.00 

1112.27 

0.00 

3.92 

0.00 

0.00 

FLOW      878. 

N  VALUES     .060 

.072  C 

.000 

RCH  = 

9300. 

10 

W      .0030261 

0.00 

3.67 

0.00 

.21 

1138.46 

1.04 

0. 

238. 

0. 

0.00 

1138.46 

0.00 

2.91 

0.00 

0.00 

FLOW      572. 

N  VALUES     .071 

.085  0 

.000 

RCH  = 

5900. 

9 

W      .0011527 

0.00 

2.50 

0.00 

.09 

1147.89 

.81 

0. 

228. 

0. 

0.00 

1147.86 

0.00 

3.38 

0.00 

0.00 

FLOW      536. 

N  VALUES     .105 

.085 

.  150 

RCH  = 

6000. 

8 

W      .0024488 

0.00 

2.12 

0.00 

.07 

1162.48 

1.06 

0. 

252. 

0. 

0.00 

1162.49 

0.00 

1.50 

0.00 

0.00 

FLOW      388. 

N  VALUES     .040 

.070 

.  150 

RCH  = 

3700. 

7 

W      .0034720 

0.00 

3.13 

0.00 

.15 

1173.83 

1.07 

0. 

123. 

0. 

0.00 

1173.84 

0.00 

2.07 

0.00 

0.00 

FLOW      361. 

N  VALUES    0.000 

.065 

.040 

RCH  = 

3700. 

6 

W      .0029596 

0.00 

2.70 

0.00 

.11 

1184.14 

.99 

0. 

133. 

o. 

0.00 

1184.  14 

0.00 

1.87 

0.00 

0.00 

FLOW      349. 

N  VALUES     .082 

.055 

.070 

RCH  = 

6900. 

5 

W      .0040078 

0.00 

2.91 

0.00 

.13 

1211.12 

1.04 

0. 

119. 

o. 

0.00 

1211.10 

0.00 

1.67 

0.00 

0.00 

FLOW      245. 

N  VALUES     .075 

.050 

.075 
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TABLE   6    (CONTINUED) 
RCH=   5100. 

4  W      .0036775 
1.01 

FLOW      232. 
N  VALUES    C.000 
RCH=   3900. 

3  W      .0041448 
1.02 

FLOW      138. 
N  VALUES    0.000 
RCH  =   6500. 


0.00 

2.31 

0.00 

.08 

1230.01 

0. 

105. 

0. 

0.00 

1230.03 

0.00 

1.26 

0.00 

0.00 

.060 


.060 


,040 


0.00 

2.57 

0.00 

.10 

1246.49 

0. 

90. 

0. 

0.00 

1246.50 

0.00 

1.35 

0.00 

0.00 

,065 


2    W               .0045985 

0.00             2.52 

0.00 

.09 

1275.74 

1.04 

0.                54. 

0. 

0.00 

1275.75 

0.00             1.21 

0.00 

0.00 

FLOW                103. 

H    VALUES             .050 

.055    0.000 

RCH=       800C 

1    W                .0047852 

0.00             2.61 

0.00 

.10 

1313.62 

.94 

0.                39. 

0. 

0.00 

1313.66 

0.00             1.24 

0.00 

0.00 
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TABLE 


STORAGE  INDICATION  MACHINE  OUTPUT 


TATION 

CFS 

S 

DELT 

(0/2  +  S/DELT) 

AREA 

22 

W 

3173. 

.38133369E+04 

.25 

.16840247E+05 

926. 

21 

W 

3107. 

.26715885E+04 

.25 

.12240004E+05 

891. 

20 

w 

3042. 

.10933095E+04 

.25 

.58944380E+04 

940. 

19 

w 

2986. 

•12969896E+04 

.25 

.66813584E+04 

573. 

18 

w 

2118. 

.11332441E+04 

.25 

.55921764E+04 

593. 

17 

w 

2092. 

.33452361E+04 

.25 

. 14427344E+05 

606. 

16 

N 

2050. 

.18138548E+04 

.25 

.82804192E+04 

1204. 

15 

W 

1930. 

•11042246E+04 

.25 

.53e20484E+04 

427. 

14 

w 

1097. 

.74287252E+03 

.25 

.35204400E+04 

336. 

13 

w 

1040. 

.44508772E+03 

.25 

.23005008E+04 

307. 

12 

w 

979. 

.10612625E+04 

.25 

.47350000E+04 

373. 

11 

w 

946. 

.69875544E+03 

.25 

.32682217E+04 

302. 

10 

w 

878. 

.38267458E+03 

.25 

.19697483E+04 

238. 

9 

w 

572. 

.4006C011E+03 

.25 

.18884504E+04 

228. 

8 

w 

536. 

.19335186E+03 

.25 

.10418074E+04 

252. 

7 

w 

388. 

.13224861E+03 

.25 

.72329444E+03 

123. 

6 

w 

361. 

.24280616E-t03 

.25 

.  11521246E+04 

133. 

5 

w 

349. 

.15980334E+03 

.25 

.81406336E+03 

119. 

4 

w 

245. 

.10636202E+03 

.25 

.54809808E+03 

105. 

3 

w 

232. 

.13131788E+03  ' 

.25 

•64172152E+03 

90. 

2 

w 

138. 

.10510562E+03 

.25 

.48977248E+03 

54. 

1 

w 

103. 

39. 
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TABLE   8       LISTING  OF  SOURCE  PROGRAM   ARS   USHL   2. 

DIMENSION  X<32)»  Y(32).  C < 3 ) . A ( 3 ) »P < 3 ) »R ( 3 ) . DELQ < 3 ) . V20 ( 3 ) . V ( 3 ) . RR 

1 

1(3) »AA(3) »CC(3) 

2 

COMMON  A.P.N.EWS.OBL.OBR.X.Y.NSTA.LET.LL.YDELO.NPASS.AN.RCH.SCH.RC 

3 

1HO.EWS2.EVRT.SLOPOO 

4 

235 

ECOMO=0. 

5 

ECOMP  =0. 

6 

YDELO  =  0. 

7 

HVO  =0. 

8 

HBR  =  0. 

9 

CHBO=0. 

10 

SLOPO=0. 

11 

K  =  1 

12 

LN  =  0 

13 

KODE=l 

14 

PAUSE 

15 

IF(SENSE  SWITCH  3)515*516 

16 

515 

PUNCH  530.KODE 

•17 

516 

LN1  =  0 

18 

READ  220,EI»CSM.SLNEW»HVNEW 

19 

NCSM=CSM 

20 

PUNCH  230»NCSM 

21 

99 

READ  l»NSTA»LET»RCH,OBL.OBR»YAR<»BEND*SQZE  »EVRT»FLAG 

22 

BENO  =  0. 

23 

SQZO  =  0. 

24 

EINT=EI 

25 

RCHO=RCH 

26 

ECOMOO=ECOMO 

27 

SLOPOO=SLOPO 

28 

AN=1. 

29 

IF(FLAGJ2»4.2 

30 

2 

READ  3. DA. (C( I ) ♦ I =1 ♦ 3 > »ELEV »DT 

31 

IF(DT)520. 520.521 

32 

521 

DELT=DT 

33 

520 

FLOW=DA*CSM 

34 

FLOW2=FLOW 

35 

IFIELEV195.150.95 

36 

95 

K    =  1 

37 

PUNCH  58.  FLOW.ELEV 

38 

154 

PUNCH  59.  (C(  I  )  .1  =  1.3) 

39 

GO  TO  4 

40 

150 

PUNCH  152. FLOW 

41 

GO  TO  154 

42 

4 

DO  6  I  =  1.29.4 

43 

11=1+3 

44 

READ  5.NSTA.LET.(X(J).Y(J).J=I»II).EX 

45 

M  =  1  +  3 

46 

DO  102  N=I .M 

47 

IF  <Y(N>)  102.7.102 

48 
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TABLE   8    (CONTINUED) 

102 

CONTINUE 

49 

101 

IF(EX)7.6»7 

50 

6 

CONTINUE 

51 

7 

N  =  N-1 

52 

CHB=Y( 1 ) 

53 

DO  130  I = 1 »N 

54 

IF(CHB-Y(  1)1130.130.120 

55 

120 

CHB=Y( I ) 

56 

130 

CONTINUE 

57 

IF<  LN)  110.205.110 

58 

205 

LN=1 

59 

GO  TO  139 

60 

110 

CHE=CHB-CHBO 

61 

IF(CHE)  111.111.112 

62 

111 

CHE=0. 

63 

112 

SCH=CHE/RCHO 

64 

209 

IFtSENSE  SWITCH  1)215.226 

65 

215 

PRINT  217 

66 

ACCEPT  220.EINT 

67 

RCH=RCHO 

68 

ECOMO=ECOMOO 

69 

SLOPE=SLOPOO 

70 

SLOPO=SLOPE 

71 

AN=1. 

72 

226 

NSECS=CHE/EINT 

73 

IF (NSECS-3 1227.22 8. 228 

74 

227 

NSECS=3 

75 

228 

CFAC=NSECS 

76 

RCH=RCH/CFAC 

77 

GO  TO  140 

78 

139 

CFAC=1 

79 

140 

GO  TO( 10.9) »K 

80 

C   TEST  FOR  BRIDGE  STATION  FLAG 

81 

9 

IFtEVRT 111.12.11 

82 

11 

EWS  =  ECOMO 

83 

GO  TO  302 

84 

12 

EWS  = ( SLOPE»RCH)+ECOMO 

85 

GO  TO  13 

86 

10 

EWS  =  ELEV 

87 

302 

EWS2=EWS 

88 

13 

DO  14  J  =  1.3 

89 

A(J)  =  0. 

90 

P( J)  =  0. 

91 

R( J)  =  0. 

92 

V(  J)  =  0. 

93 

DELQfJ)  =  0. 

94 

14 

V2GU)  =  0. 

95 

LL  =  0 

96 

CALL  AREA 

97 

IF(LL-1  129.99.209 

98 

29 

IF(EVRT)  30.32.30 

99 

30 

OMEGA  =  HVO  /(EWS  -  EVRT) 

100 

ALFA  =  (ASUM  -  A(l)  -  A(2)  -  A(3))  /  ASUM 

101 

30 


TABLE   8    (CONTINUED) 

AL4  =  ALFA  »  ALFA  *  ALFA  *  ALFA  102 

HBR  =  2.*  YARK*(YARK  +10 .*OMEGA- .6 ) * ( ALFA+ 15. *AL4 )  *HVO  103 

PUNCH  31 .NSTA.LET. ( A(  I  )  .1=1.3)  »HBR  104 

GO  TO  99  1°5 

32  QSUM=0.  106 
DO  34  J=1.3  107 
IF(P( J) )605.605.33  108 

33  R( J)=A( J)/P( J)  1°9 
605  IF( LN1) 182.182.180  11° 

180  A( J)=(A( J)»AN+AA< J)»(CFAC-AN) ) /CFAC  111 
R( J)=(R( J)*AN+RR< J)*(CFAC-AN) ) /CFAC  112 
C(J)=(C(J)*AN+CC( J)*(CFAC-AN) ) /CFAC  113 
IF(J-2) 182.600.182  114 

600  FLOW=( FLOW2*AN+FLOWl*(CFAC-AN) )/CFAC  115 

182  IF(C( J)  )34.34.181  116 

181  V( J)=.01486*(R( J)**.67)/C( J)  117 
DELQU)  =  V(J>  *  AtJ)  118 

34  QSUM  =  DELQ(J)  +  QSUM  119 
RATIO  =  FLOW/  QSUM  120 
SLOPE  =  RATIO  *  RATIO  *  .0001  121 
GO  TO  (36.35) »K  122 

35  HFRIC  =  (SLOPE  +  SLOPO)  »  .5  »  RCH  123 
GO  TO  37  124 

36  HFRIC  =  0.  125 

37  SUMVQ  =  0.  126 
DO  39  J  =  1.3  127 
IF  (V( J)  J96.96.38  128 

96  V(J)  =  0.  129 
V2Q(J>  =  0.  130 
GO  TO  39  131 

38  V(J)  =  V(J)  *  RATIO  I32 
V20U)  =  V(J)  *  V(J)  *  DELQ(J)  *  RATIO                                    133 

39  SUMVQ  =  SUMVQ  +  V2Q(J)  13A 
IF(HVNEW)621.621.620  I35 

620  HV1=HVNEW  I36 
HVNEW=0.  I37 
GO  TO  625  I38 

621  HV1=SUMVQ/(FL0W*64.32)  I39 
625  GO  TO(97.100) .K  1^0 

97  ECOMP  =  ELEV  I41 
GO  TO  49  I42 

100  DELHV  =  HV1  -  HVO  l^3 

IF(DELHV)41. 44.40  1** 

41     DELHV    =    -    DELHV  1^5 

40  IF(SQZE )44. 44.700  1^6 

700  IF(  AN-CFAO44.80.80  1*7 
80  SQZO  =  DELHV  *  SQZE  1*8 
44     IF(BEND)60. 60.701  1*9 

701  IF(  AN-CFAO60.82.82  150 
82  BENO  =  (HV1  +  HVO)  *  .5  *  BEND  151 
60    ECOMP    =    ECOMO    +    HFRIC    +    SQZO    +  BENO    +    HVO    -    HV1                                                                    152 

ETEST  =  (ECOMP  -  EWS)  -.1  153 

IF(ETESJ)  47.  49.48  154 
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TABLE   8    (CONTINUED) 

47 

ETEST  =  EWS  -  ECOMP  -.1 

155 

IF(ETEST)  49»  49.  48 

156 

48 

EWS  =  ECOMP  +  (EWS  -  ECOMP)  *  .5 

157 

GO  TO  13 

158 

49 

IF(SLNEW)301. 301.300 

159 

300 

SLOPE=SLNEW 

160 

SLNEW=0. 

161 

301 

SL0P0=SL0PE 

162 

HVO  =  HV1 

163 

ECOMO  =  ECOMP 

164 

CHB0=CHB 

165 

ASUM  =  All)  +  A(2)  +  A(3) 

166 

IF(HBR) 57.57.81 

167 

81 

IF(  AN-CFA083.55.55 

168 

55 

ELEV  =  ECOMP  +  HBR 

169 

HBR=0. 

170 

K  a  1 

171 

GO  TO  10 

172 

57 

HBR=0. 

173 

83 

K  =  2 

174 

AN=AN+1. 

175 

NPASS=0 

176 

IF(AN-CFAC)9.650.28  0 

177 

650 

IF(EWS-CHB)651.9.9 

178 

651 

EWS=CHB 

179 

GO  TO  13 

180 

280 

IF( LN1 1499.499.500 

181 

499 

LN1  =  1 

182 

GO  TO  510 

183 

500 

IF(SENSE  SWITCH  3)505.510 

184 

505 

A1  =  0. 

185 

A2  =  0. 

186 

DO  507  1=1.3 

187 

A1=A1+A( I ) 

188 

507 

A2=A2+AA( I ) 

189 

AVGA=( A1+A2J/2. 

190 

S=RCHO»AVGA/3600. 

191 

R0F1  =  ( FLOW  1/2. 1+S/DELT 

192 

PUNCH  5  40.NSTAO.LET.FLOW1 »S . DELT .ROF 1 . A2 . KODE 

193 

IF(NSTA-1 ) 510.508.510 

194 

508 

PUNCH  541 .NSTA.LET.FL0W2.A1.K0DE 

195 

510 

DO  85  J=l .3 

196 

RR( J)=R( J) 

197 

AA( J)=A( J) 

198 

85 

CC( J  )=C( J) 

199 

FL0W1=FL0W2 

200 

NSTAO=NSTA 

201 

RCH=RCHO 

202 

PUNCH  52. RCH.NSTA .LET. SLOPE. <V<I).I=1.3) »H VI. ECOMP 

203 

PUNCH  53,HFRIC.(A(I).I=1.3).BENO.EWS.(R(I).I=1.3) »SQZO 

204 

IF(NSTA-1 )99.240.99 

205 

240 

PRINT  234.NCSM 

206 

PAUSE 

207 
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TABLE   8    (CONTINUED) 

GO  TO  235  208 

1  FORMAT! I4»A2.F5.0*2F4.0»3F3»2»F5.1»F1.0)  209 

3  FORMAT!6X»F6.2»3F3.3.F6.2»F4.2)  210 

5  FORMAT!  I 4» A2 > 8F5. 1 » Fl .0 >  211 

31  FORMAT! 14HBRIDGE  SECT  ION / I  5 ♦ A2 ♦ 1  OX , F 10,0 » 2F8 . 0 » F 10 . 2  )  212 

58  FORMAT!6X4HFLOW.F9.0»5X19HSTARTING  W.S.  ELEV..F9.2)  213 

52  FORMAT(4HRCH=,F7.0/I5»A2»4X»F9.7  »1X,3(2X»F6.2)»3X,F5.?»3X»F8.2>  214 

53  FORMAT(9X.F6.2»2X»F10.0»2F8.0»F10.2,F11.2/21X.4F8.2 )  215 

59  FORMAT! 6X»8HN  VALUES  * 2X » 3F6 . 3 )  216 
152  FORMAT(6X4HFLOW»F9.0)  217 
217  FORMAT! 23HENTER  NEW  ELEV  INTERVAL)  218 
220  FORMAT! F 3 . 2 » F6 . 0 » F7 . 7 » F4 . 2)  219 
230  FORMAT! I5»12H  CSM  PROFILE//)  220 
234  F0RMATI6HEND  OF.I5.49H  CSM  PROFILE-LOAD  DATA-PUSH  START  FOR  NEXT  P  221 

IROFILtIHE)  222 

530  F0RMATI49H  STATION    CFS  5  DFLT  (0/2  +  S/»14HDEL  223 

IT)       AREA.15X.I2)  224 

540  FORMAT!  I  5 ♦ A2 »F9 . 0 » 1 X » E 14 . 8 » F8 . 2 ♦ IX » E 1 4 . 8 » F  1 1 . 0 » 1 3X » 12 )  225 

541  FORMAT!  I  5 » A2 ♦ F9. 0 ♦ 38X . F 1 1 . 0 ♦ 1 3X ♦ 12)  226 
END  227 


SUBROUTINE  AREA  1 

DIMENSION  A! 3) »P< 3) »X( 32) tY<32)  2 

COMMON  A,P,N»EWS»OBL.OBR»X,Y,NSTA»LET.LL.YDELO,NPASS»AN,RCH,SCH,RC  3 

1H0.EWS2»EVRT»SL0P00  4 

XBASF(E»G)=E-G  5 

XYARF!E»G»H)=E*G/H  6 

XATSF(E»G»H)=E*<G-H)  ? 

XWPPF(E*G»H,D)=SQRTF! ( E-G ) **2+ < H-D ) **2 )  8 

IF(RCHO)205.205.197  9 

197  IF(EVRT)205.199,205  10 

199  EWS2=(RCHO-AN*RCH)*SCH+EWS  11 
IF(EWS2-Y( 1)  )200.200.201  12 

200  IF(EWS2-Y(N) )205»205»202  13 

201  LM=1  14 
GO  TO  203  15 

202  LM=N  16 

203  NPASS=NPASS+1  I7 
IF! NPASS-2)204»280»280  18 

204  EWS2=Y(LM)  19 

205  YDELO=0.  20 
DO  270  1  =  1. N  21 
IF(OBL-X(  I  )  )212»210»210  22 

210  J=l  23 

GO  TO  220  24 
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TABLE   8    (CONTINUED) 

212  IF(OBR-X<  I  )  )216»215.215  25 

215  J=2  26 
GO  TO  220  27 

216  J=3  28 
220  YDEL  =  XBASF(EWS2»Y(  I  )  )  29 

IF( YDEL)225.225t230  30 

225  IF< YDELO1270. 270.250  31 

230  IF(YDELO)235. 265.258  32 

235  K=l  33 

236  B=XBASF(X<  I )  »X<  1-1 )  )  34 
C=XBASF( Y( I) »Y( 1-1)1  35 

2^0  GO  TO(245. 255.260. 264.267) .K  36 

245  XWS=X( I l+XYARFI B.YDEL.C)  37 

A< J)=A< J)+XATSF( YDEL/2.»X(  I  )  ,XWS  )  38 

P( J) =P( J)+XWPPF(X( I ) , XWS.YDEL.O.  )  3  9 

GO  TO  270  40 

250  K=2  41 

GO  TO  236  42 

255  XWS=X( I )+XYARF(B.YDEL»C>  43 

A(J)=A(J)+XATSF<YDELO/2.»XWS»X(I-l))  44 

P( J) =P( J)+XWPPF(XWS.X (I-l).YDELO.O.)  45 

GO  TO  270  46 

258  IF(Y( I )-Y( 1-1 )) 259.263.263  47 

259  K=3  48 
GO  TO  236  49 

260  A( J)=A( J)+XYARF(B.C.-2. ) +XYARF < B ♦ YDELO . 1. )  50 
P( J)  =P( J)+XWPPF(B.O..C.O.  )  51 
GO  TO  270  52 

263  K=4  53 
GO  TO  236  54 

264  A( J) =A( J1+XYARFIB.C.2. ) +XYARF < B . YDEL » 1 . )  55 
P< J) =P( J)+XWPPF(8.0..C.O.  )  56 
GO  TO  270  57 

265  K=5  58 
GO  TO  236  59 

267  At  J) =A( J1+XYARFIB.C.-2. )  60 

P( J)=P( J)+XWPPF(B.O..C.O.  )  61 

270  YDELO=YDEL  62 

GO  TO  285  63 

280  NPASS=0  64 

PRINT  300.NSTA.LET  65 

PAUSE  66 

IFISENSE  SWITCH  1)283.284  67 

283  LL  =  2  68 
GO  TO  285  69 

284  LL=1  70 

285  RETURN  71 
300  FORMATI21HEXCEEDS  SURVEY  L I M I TS . 1 4 . A2 » 2X .49HTURN  SW  1  ON  FOR  RERUN  72 

1  WITH  SHORTER  REACH  LENGTH  /33HOFF  CONTINUES  UPSTREAM .PUSH  START)  73 

END  T* 
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JDe-tailed   Flow  Diagram- -Computer   Program   USHL  #2 


/fcEAD.EI  .CSH,   "\ 

VsLSEV , HVNEW      J 


c 


)ncsm=csm1 

Pl^NCH  ,NCSM 


J 


f~~\       ,  /read, s?ta, ut, rch.obl.obr^ 

(    5     )         MyaR*  , BESD .  SQZE ,  EVRT ,  FLACy 


BENO^O . 
SQZO=0 . 
EINT=EI 
RCHO-RCH 
EC0M0O=ECOMO 
SLOPOO=SLOPC 
SLOPO   "   SLCPE 
AN   -   1. 


CREAD,  NSTA,  LET,  (X(J) ,       *\ 
Y(J),J=I,II),EX  J 


~\        /ACCEPT,  EIlfT     \ 


RCH=RCHO 
ECOMO-ECOMOO 
SLOPE=SLOPOO 
A.1-1. 
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Detailed   Flow  Diagram   (continued) 


OMEGA=HVO/  (EWS-EVKT) 
ALFA»(ASUM-A(1)-A(2)-A(3))/ASUM 
ALA=ALFAxALFAxALFAxALFA 
HBR=2.x*AR*x(YARK+10.xOMECA-.6)x(ALFA+15.xALA)xHV0 


(PUNCH,NSTA,LBT,        ^^         /p 
(*(I),I»1,3),HBR  J  \^_J 


A(J)=  (A(J)xAH+AA(J)x(CFAC-AlO)  /CFAC 
R(J><R(J)xAN+RR(J)x(CPAC-AN))/CFAC 
C(J)=(C(J)xAN+CC(J)x(CPAC-AN))/CFAC 


| BENP°(HVl+HVO)x.SxBEND 


ECOMP^ECOKO+flFRIC+SQZO+BENO+flVO-HVl 
ETES  C= (ECOMP-EWS) - . 1 


£ 


ETEST=EWS-ECOKP-.l 


EWS=ECOMP+(EWS-ECOMP)x. 5 


I  HFRIP(SL0?E+SU>PO)x.5xRCB  | 


X 


1 


SOMV(^0. 


SLOPE=SLNEW 

SLNEW=0. 


SL0PO=SL0PE 

HVO=HVl 

ECOMCbECOMP 

CHBO=CHB 

ASUK=A(1)+A(2)+A(3) 
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Detailed   Flow  Diagram   (continued) 


AVGA=(Al+A2)/2. 
S=RCHOxAVGA/  3  600 . 
R0Pl=(FL0Wl/2.)+S/DELT 


1 

C      PUNCH, NSTAO, LET, FL0Wl\ 
S,DELT,ROFl,A2,KODE    J 


c 


PUNCH,  NSTA,  LET ,  FLOW2 
A1,K0DE 


0 


t 


RR(J)=R(J) 
AA(J)=A(J) 
CC(J)=C(J) 


FLOWl=FLOW2 
NSTAO=NSTA 
RCH=RCH0 


/PUNCH,  RCH,  NSTA,  LET,  SLOPE/\ 
I      (V(I),I=1,3),HV1,EC0MP 


/PUNCH, HPRICj (A(I), 1=1, 3) , 
1     BEN0,EWS,(R(I),I=1,3) 
\ SQZO 


/'PRINT,  END  OF,NCSM,CSM     \ 
[     PROFILE-LOAD  DATA-PUSH 
I     START  POR  NEXT   PROFILE    J 
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Arithmetic  Function  Statements 

XBASF(E,G)=E-G 

XYARF(E,G,H)=ExG/H 

XATSF(E,G,H)=Ex(G-H) 

XVPPF(E,G,H,D)=SQRIF((E-G)xx2+(H-D)xx2) 
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Detailed  Flow  Diagram  (continued) 

® 


XWS.X(I)+XYARF(B,YDEL,C) 

A(J)-A(J)+XAISF<Yr>EL/2.,X(I),XWS) 

P(J)=.P(J)+XVPPP(X(I),XWS,YDEL,0.) 


XyS=X(I)+XYARP(B,TOEL,C) 

A(J)-A(J)+XAISF(YDEL0/2.,XWS,X(I-1)) 

P<J)»P(J)+XVPPP(XWS.X(I-1),YDELO,0.) 


<$_ 


PRINT, EXCBEDS   SURVEY 
LIMTTS,NSTA,  LET, TURN 
SV   I   ON   POP.   RERUN  WITH 
SHORTER   REACH  LENGTH 
OFF   CONTINUES   UPSTREAM, 
PUSH   START 


©h: 


-A   LUi      h— (») 


A(J)  =  A(J)+XYARP(B,C,-2.) 
P(J)=P(J)+XWPPF(B,0.,C,0.) 


A(J)"A(J)+XKAR7(B,C,2.)+XKARf'(B,YDEL,l.) 
P(J)  =  PU)+XVPPF(»,0.  ,C,0.) 


A(J).A(J)+XYARF(B,C,-2.)+XYARP(B,YDELO,l.) 
P(J)-P(J)+XVPPP(B,O.,C,0.) 
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